Precision measurement of the quasi— free pn-^-d<fi reaction close to threshold 
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The quasi-free pn—>d(j> reaction has been studied at the Cooler Synchrotron COSY- Jiilich, using 
the internal proton beam incident on a deuterium cluster-jet target and detecting a fast deuteron 
in coincidence with the K + K~ decay of the 0-meson. The energy dependence of the total and 
differential cross sections are extracted for excess energies up to 80 MeV by determining the Fermi 
momentum of the target neutron on an event-by-event basis. Though these cross sections are 
consistent with s-wave production, the kaon angular distributions show the presence of p waves at 
quite low energy. Production on the neutron is found to be stronger than on the proton but not by 
as much as for the 77-meson. 



PACS numbers: 25.40.Ve, 13.75.Cs, 14.40.Cs 

Meson production provides access to the internal struc- 
ture of baryons and the dynamics of hadronic reac- 
tions and thus is an important exploration field for non- 
perturbative QCD. In proton-proton collisions, meson 
production has been extensively studied and data are 
now available on the production of most members of the 
fundamental pseudoscalar and vector nonets near their 
respective threshold 0, including the tt(140), 77(547), and 
77' (958), as well as the p(770), a; (782), and 0(1020). The 
</>-meson is of particular interest because of its compar- 
atively large mass and its dominant ss quark structure. 
However, in order to study all facets of meson produc- 
tion dynamics, it is necessary to investigate the isospin 
dependence by precision measurements in both pn as well 
as pp collisions. In the case of the 77-meson, such experi- 
ments have revealed that the pn production cross section 
is over six times larger than that for pp . Analogous <f> 
data are important for nucleon-nucleon production mod- 
els and also serve as crucial input in the interpretation of 
nucleon-nucleus, and nucleus-nucleus results, where in- 
medium effects are anticipated 0. 

In the absence of a free neutron target or a qual- 
ity neutron beam, quasi-free production on deuterium 
pd^dXp S p has often been substituted. Here the re- 
action is assumed to have taken place on the neutron 
bound in the deuteron and p sp is a slow "spectator" pro- 
ton that does not take an active part in the reaction 
and whose momentum reflects the Fermi motion of the 



particle before the production. In order to show that the 
reaction involved only the neutron, the spectator must be 
identified and the precise determination of the cm. en- 
ergy requires that the p sp momentum is well measured. 
Spectators emerging from an ultra-thin target with a few 
MeV can be studied directly at a storage ring using solid- 
state counters, as has been done for pd — > dir° p sp [4| and 
pd^ duj p sp ,5] . The alternative approach is to identify 
the produced meson X through its decay products and 
then reconstruct the spectator momentum using kine- 
matics. This method has been successfully employed for 
the pd—>di]p sp and pd~^pnrjp S p reactions, where the n 
was identified through its 27 decay branch ,2j. We have 
studied for the first time quasi-free pn — ► d<p production 
through the indirect method of measuring the spectator 
momentum using the K + K~ decay of the <p in coinci- 
dence with a fast deuteron. 

The experiment was performed with a 2.65 GeV pro- 
ton beam at an internal target station of the Cooler Syn- 
chrotron COSY, employing the magnetic spectrometer 
ANKE 6] to identify and measure the reaction. ANKE 
has detection systems placed to the right and left of 
the emerging beam to register slow positively and nega- 
tively charged ejectiles, with fast positively charged par- 
ticles being measured in the forward system. The deu- 
terium cluster-jet target Q provided areal densities of 
~ 3.4 x 10 14 cm~ 2 which, combined with a typical proton 
beam intensity of ~ 6.2 x 10 16 s _1 , gave an integrated 
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luminosity of 23 pb" 1 over the 300 hours of data taking. 

The pd — > d(j)p sp reaction was studied in a manner anal- 
ogous to that successfully employed for the pp — > pp<f) 
reaction at COSY ||, using the — > K + K~ decay. 
Charged kaon pairs were detected in coincidence with a 
forward-going deuteron, requiring that the overall miss- 
ing mass in the reaction was consistent with that of the 
non-observed slow spectator proton p sp . As a first step, 
positive kaons are selected through a procedure described 
in detail in Ref. |2|, using the time of flight (TOF) be- 
tween START and STOP scintillation counters of a dedi- 
cated K + detection system. In the second stage, both the 
coincident K~ and forward-going deuteron are identified 
from the time-of-flight differences between the STOP 
counters in the negative and forward detector systems 
with respect to the STOP counter in the positive system 
that was hit by the K + . These two TOF selections, as 
well as that for the K + , were carried out within ±3er 
bands. 

Fig. which shows the missing mass spectrum as- 
suming that the detected particles are indeed K + , K~ 
and deuteron, demonstrates a clear peak at the mass of 
the proton. The secondary peak around 1.02 GeV/c 2 
arises from pir + Tr~ events, where a n + was misidenti- 
fied as a K + . This background is well separated from 
the spectator peak over the whole kinematic region. The 
residue from misidentified particles inside the proton gate 
(of ± 3 a) is 3.1 % and such events generally fail the later 
criteria of the analysis. In total, about 4500 coincidences 
were retained as dK + K~p sp events for further study. 

The K + K~ invariant -mass spectrum for the 4500 
events is shown in Fig. \\]p. The distribution is domi- 
nated by the </()-meson peak, which sits on a slowly vary- 
ing physical background from direct K + K~ production. 
This has been estimated by a four-body phase-space sim- 
ulation which, together with the 4> contribution, is fitted 
to the overall spectrum. The shape of the resonant con- 
tribution is reproduced by the natural width of 0-meson 
with an experimental mass resolution a = 1 MeV/c 2 , 
which is consistent with the momentum resolution of 
the ANKE detector system. The direct K + K~ con- 
tribution, which is less than 8% in the </> mass region 
1.020 ± 0.015 GeV/c 2 , could be easily subtracted. 

The momentum distribution of the unobserved proton 
for events in the <fi peak is shown in Fig.^Ji. As expected 
for a spectator proton, this spectrum is peaked at very 
low values and there are few events with momenta above 
about 150 MeV/c. To confirm the spectator hypothesis, 
a Monte Carlo simulation has been performed where the 
Fermi momentum in the target deuteron has been derived 
from the Bonn potential 10] . The energy dependence of 
the pn — ► d(p cross section is assumed to follow phase 
space, which is consistent with the results to be shown 
later. After including the detector response, the simu- 
lation fits very well the shape of the data for momenta 
up to at least 150 MeV/c, a region where the model de- 
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FIG. 1: a) Missing-mass distribution of the pd— > dK + K~ X 
events with lines indicating the proton selection range, b) 
K + K~ invariant -mass distribution with lines showing the 
mass range of selected 0-meson events. The dashed curve 
shows the four-body phase-space fit to the non-resonant pro- 
duction while the solid histogram is the sum of this and the <j>- 
meson contribution, c) Momentum distribution of the unob- 
served proton compared with a Monte-Carlo simulation based 
on the spectator model, d) Distribution in cm. excess energy 
compared with the same simulation. 



pendence of the deuteron wave function is negligible com- 
pared with our statistical uncertainty. The spectator dis- 
tribution could be obtained with even greater precision 
than that for the corresponding pd — > drjp sp reaction |2|. 

Due to the Fermi motion of the neutron in the 
deuteron, the cm. excess energy e = ^fs — (m c i + m<p) 
is spread over a range of values even for a fixed beam 
energy. Since we have completely determined the kine- 
matics for each of the pd—>d(f>p sp events, the value of e 
could be calculated on an event-by-event basis, and the 
resulting distribution is shown in Fig.^i. This is also well 
described by the simulation, which shows that e can be 
reconstructed with an average precision of a e = 2 MeV, 
and which can be used in the extraction of cross sections 
for e < 80 MeV. 

The target density was determined by measuring the 
frequency shift of the stored proton beam as it lost en- 
ergy due to its repeated passages through the target 01 • 
Combined with measurements of the beam current this 
yielded the value of the luminosity L with a precision of 
about ±6 %. This was checked through the simultaneous 
measurement of pd elastic and quasi-elastic scattering, 
where a fast proton was registered in the polar angular 
range 5.0° < ^ < 8.5° in the forward detectors. The lu- 
minosity was then obtained from estimates of the cor- 
responding differential cross sections within the Glauber 
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FIG. 2: Differential cross sections for pn — > d(j> for different 
ranges of excess energy, with the left panels showing the de- 
pendence on the angle of the K + from the 4> decay in the (ft 
rest frame, and the right panels the dependence on the polar 
angle of the <j> in the overall cm. system. Vertical error bars 
indicate purely statistical uncertainties whereas the horizon- 
tal ranges reflect the bin width. 

formalism 0] . This technique has been used successfully 
at other energies ^4| . Though the two methods give con- 
sistent results to within 3%, the error in the pd technique 
is about 10%, due mainly to the use of the theoretical 
model and uncertainty in the acceptance correction. 

In order to evaluate differential cross sections, the 
geometrical acceptance, resolution, detector efficiency 
and kaon decay probability were taken into account in 
a Monte Carlo simulation, using the GEANT4 pack- 
age (l5| . For a given excess energy, the distributions in 
all variables were consistent with phase space except for 
that of , which is the polar angle of the K + from the 
4> decay with respect to the beam direction, in the <f> rest 
frame. At threshold the only allowed pn — > d<f> transition 
arises from an initial 1 P\ state. The unique production 
amplitude is therefore of the form M — p- (e\ x et) $ pn , 
where p is the beam momentum, Ed and are the po- 
larization vectors of the deuteron and <j> respectively, and 
&pn represents the spin-0 initial pn state. From the 
structure of the matrix element, it follows that the 4>- 
meson is aligned transversally to the beam so that, fol- 
lowing its decay, the kaons cannot be produced along the 
beam direction and a sin 2 behavior is to be expected. 

To allow for the possibility of higher partial waves, the 
distribution was parameterized in the most general al- 
lowed form: dcr/dfi^ = 3(asin 2 9^ + 2fecos 2 9^)/87r, 
normalized such that the total cross section a = a + b. 
This form was handled iteratively in the simulation to 
get the best values of the parameters a and b and of the 
ANKE acceptance. For large e the acceptance in the 
backward cm. hemispheres are somewhat higher than in 
the forward, but all distributions are completely consis- 
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FIG. 3: Total cross section for the quasi-free pn — > d<f> re- 
action as a function of the excess energy (filled circles). In 
addition to the statistical error bars shown, there are over- 
all systematic uncertainties of ±10%. The curve represents a 
phase-space yfe behavior. For comparison, we also show as 
open circles the values obtained from pp — > ppcj> [g ■ 

tent with them being symmetric in the cm. system. The 
results for different excess energy bins in Fig. [21 are there- 
fore shown as functions of the magnitudes of the cosines 
of and the polar angle Of m of the (j), for which the 
resolutions are estimated to be 0.024 and typically 0.02- 
0.04 respectively The dominance of the sin 2 term is 
very clear at the lower energies and all the data are well 
represented by b/a w (0.012±0.001) (e/MeV). Given the 
ambiguity associated with the nine possible p-wave am- 
plitudes, this ratio represents the minimum fraction of 
higher partial waves and indicates that this is significant 
for the larger e. Despite this, the angular distribution of 
<f> production in the overall cm. system is consistent with 
isotropy for all e. Note that the production of p-wave <f> 
from an initial 3 Si state would also be flat in cosQf m . 

The values of the a and b coefficients lead directly 
to the total cross section for <f> production shown in 
Fig. |31 while numerical values will be found in the 
HEP database [lif. In addition to the point-to-point 
statistical errors, there is an overall systematic uncer- 
tainty of ±10% coming from luminosity (6%), stability of 
the data-taking efficiency (4%), background (3%), and 
MWPC efficiency corrections for kaon detection (4%). 
The results have not been corrected for the reduction in 
the incident flux due to shadowing by the proton in the 
deuteron target, which would increase the cross sections 
by about 4% pjj. Values of the pp — ► pp4> total cross 
sections available in our energy range are also shown || . 

Two-body phase space increases like y/e and this is 
distorted by less than 4% when the width of the 4>- 
meson is taken into account. As shown in Fig. we 
find that a{pn -> d<j>) = (48 ± 1) y/(e/MeV) nb, de- 
spite the decay angular distributions showing significant 
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p-wave effects at higher e. The values are much higher 
than those of pp — > ppcp 0, E3 , but this is due in part 
to there being a three- rather than a two-body final 
state. However, very near threshold, isoscalar S'-wave 
4>pn production can be estimated from our dip data using 
final-state-interaction theory 19| , a technique that has 
been tested for n production 2j. This approach yields 
a(pn — ► pn<p)/a(pp — » pp<p) ~ 2.3 ± 0.5, which is only 
about a third as big as the ratio for n production 

The ratio / u of the production of the light isoscalar 
vector mesons <p and u> in various nuclear reactions in- 
volving non-strange particles provides valuable tests of 
the Okubo-Zweig-Iizuka rule [13] ■ This rule suggests 
that, due to small deviations from ideal mixing of these 
mesons, one should have R$/ u « Rozi = 4.2 x 10 -3 j2l| 
under similar kinematic conditions. Significant enhance- 
ments of this ratio have, however, been reported in the 
literature and for proton-proton collisions we recently ob- 
tained — > pp4>)/a(pp — ► ppco) = (3.3 ± 0.6) x 10~ 2 ~ 
8 x i?ozi M- There is a measurement of u> production 
in proton-neutron collisions at 57T^5 MeV 0] and, com- 
paring this with our data, we find that at this energy 
a{pn -> d(f>)/a(pn -> duj) = (4.0±1.9) x 10~ 2 w 9 x i?ozi- 
Though similar to the pp result, the error bar is large. 

In near-threshold production reactions, the relevant 
degrees of freedom seem to be mesons and baryons rather 
than quarks and gluons, and the predictions for <p pro- 
duction in pn — > dip are very sensitive to meson ex- 
change and nucleonic currents [2J, |23, |24( • Nevertheless, 
all three calculations yield broadly similar values for the 
pn — > d(j> total cross section, being in the 0.1-0.5 ub range 
at e w 50 MeV compared with the « 0.34 jih shown by 
our data in Fig. EI However, whereas one calculation sug- 
gests that the <p production is maximal in the forward 
direction , another predicts it to be much flatter • 
Our data in Fig.|2are consistent with isotropy. Although 
no calculations appear to exist in the literature for the 
polarization of the 0-meson in the pn — ► dip reaction, es- 
timates have been made for that of the deuteron and of 
the initial pn spin-correlation [24| ■ Both of these observ- 
ables are sensitive to effects from p-wave <p production 
but neither shows as big effects with energy as we have 
seen from the <p alignment in Fig. [21 

In summary, we have presented the first measurements 
of <p production in pn collisions that will provide im- 
portant constraints on the theoretical modeling of such 
processes. The <p polarization, as measured through its 
K + K ~ decay, shows the early onset of p-waves which are 
not apparent in the energy variation of the total cross sec- 
tion. This behavior can also not be seen in the cm. angu- 
lar distributions and this suggests that p-waves might be 
more important than previously thought in other near- 
threshold meson production. The production is stronger 



in pn collisions than in pp, though the factor is not as 
large as for r\ production. 

These data are, of course, valuable in the interpre- 
tation of 0-meson production in the collision of heavy 
ions. Further testing of the OZI rule in pn collisions will 
have to await better data on u> production. This should 
be possible in the future at the WASA at COSY facil- 
ity |25j, where photons from the to — > 7r°7 decay can be 
detected in coincidence with fast deuterons. This decay 
mode would also allow the polarization of the to to be an- 
alyzed to see if this also shows a rapid onset of p waves. 
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